Itinerant Antiferromagnetism in Infinite Dimensional Kondo Lattice 
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Highly accurate numerical results for single-particle spectrum and order parameter are obtained 
for the magnetically ordered Kondo lattice by means of the dynamical mean-field theory combined 
with the continuous-time quantum Monte Carlo method. Hybridized energy bands involving local 
spins are identified in the Neel state as a hallmark of itinerant antiferromagnetism. At the boundary 
of the reduced Brillouin zone, the two-fold degeneracy remains in spite of the doubled unit cell. This 
degeneracy results if the molecular field felt by localized spins has identical magnitude and reversed 
direction with that of conduction electrons. The persistent Kondo effect is responsible for the 
behavior. The antiferromagnetic quantum transition occurs inside the itinerant regime, and does 
not accompany the itinerant-localized transition. 



Distinction between itinerant and localized characters 
of strongly correlated electrons has been one of the most 
fundamental issues in condensed-matter physics. The 
Kondo effect plays a central role in this problem be- 
cause even a localized spin may acquire itinerant char- 
acter by coupling with conduction electrons, and may 
form a Fermi liquid. Some of recent experiment sug- 
gest that quantum phase transition between magnetic 
and paramagnetic ground states accompanies a change 
between localized and itinerant character of electrons, 
with non-Fermi liquid behavior in the vicinity of the 
transition Another experiment using the de Haas- 
van Alphen effect has probed the change of the Fermi 
surface as the external pressure drives such systems as 
Celn3 and CeRhlns across the magnetic transition 
On the other hand, recent photoemission experiment for 
CeRu2Si2 and CcRu2Si2- a; Ge, E indicates that the Fermi 
surfaces of both are essentially the same, which involve / 
electrons and are referred to as the large Fermi surface. 
The experiment is performed at temperatures above the 
Neel transition of CeRu2Si2- a; Ge : i; [2(, while CeRu2Si2 
remains paramagnetic down at least to 0.1 K. Thus the 
antiferromagnetism in CeRu2Si2- 2; Ge 2 ; seems to be itin- 
erant. Consistent understanding of this variety of phe- 
nomena is still lacking, and accurate theoretical analysis 
is necessary to understand how the magnetic order is re- 
lated to the change from itinerant to localized characters 
of electrons. 

In this paper we show that the antiferromagnetism in 
the Kondo lattice model (KLM) occurs within the itin- 
erant regime. The KLM is the simplest system that is 
capable of describing both itinerant and localized char- 
acters of electrons, and is given by 



H = ^ e k cl a c ka + 2 J Si ■ s c 



(1) 
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where the first term represents the kinetic energy of con- 
duction electrons, Si is the localized spin at site i, and 
s c i denotes the conduction-electron spin at the same site. 
In terms of creation c\ a and annihilation Ci a operators 
of conduction electrons at site i with spin a, we obtain 



IE 



a/3 ' 



r af3Ct(3- 



The KLM has been investi- 



gated beyond the mean field theory in one- 0, [B| and 
two- [1, 0, 0, @, E3] dimensional systems. We approach 
the KLM from infinite dimensions using the dynamical 
mean-field theory (DMFT) to allow for the Neel state at 
finite temperatures. We use the continuous-time quan- 
tum Monte Carlo method (CT-QMC) [U [H, [H as the 
impurity solver, and focus on the case of unit number 
n c = 1 of conduction electrons per site. In contast with 
the metallic antiferromagnetism with n c ^ 1, where an 
incommensurate order may exist, we can safely assume 
the simple staggered order because of the nesting condi- 
tion for the conduction band in the hypercubic lattice. 
Hence the results in this paper are exact in infinite di- 
mensions except for statistical errors. 

As the half-filled limit is approached from n c ^ 1, the 
large Fermi surface tends to the zone boundary of the 
paramagnetic phase, while the small Fermi surface in- 
volves half of the Brillouin zone volume. Provided that 
the ground state has no discontinuity in the zero-doping 
limit, the limiting location of the Fermi surface should 
be reflected in the location of the energy gap in the half- 
filled case. In the localized antiferromagnetism, the gap 
opens at the boundary of the new Brillouin zone since 
conduction electrons feel the staggered internal field. In 
itinerant magnetism, on the other hand, we show in this 
paper that the energy gap is located in the center of the 
new Brillouin zone. The zone-center location is due to 
emergence of energy bands of magnetic electrons. Hence 
location of the energy gap distinguishes between itinerant 
and localized behaviors. Note that such distinction does 
not apply to the single band model such as the Hubbard 
model where the energy gap in the half-filled limit always 
occurs in the boundary of the Brillouin zone. In this case 
the character of electrons changes continuously from the 
itinerant limit to the localized one as the Coulomb repul- 
sion increases relative to the band width. 
We take the bare density of states 



p c (u>) = y / 2/irexp(-2uj 2 ) 



(2) 



which corresponds to infinite-dimensional hypercubic lat- 
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tice. We have taken the band width D = 1 as the unit 
of energy. For the Gaussian density of states, the Kondo 
temperature Tk is defined by 



e -7/2 
p c (0)T K = — ^exp 



1 



2p c (0)J 



(3) 



where e~ 7/2 /V7r ~ 0.42 with 7 ~ 0.577 being the Euler 
constant. This expression of Tk corresponds to diver- 
gence of effective exchange in the lowest-order scaling 
equation. 

The hypercubic lattice has the nesting property with 
the wave vector Q = (tt, tt, ■ ■ ■ ) at half filling, which fa- 
vors the staggered order. In the two-sublattice formal- 
ism, the Green function Gka(z) of conduction electrons 
is a 2 x 2 matrix where z is a complex energy, and a 
wave vector k belongs to the reduced Brillouin zone 14 1 . 
In the DMFT, the wave vector enters only through Ek- 
Therefore we introduce the notation n = and regard 
k as if it represents a wave number. The spectral func- 
tion A a (K,uj) can be calculated from the matrix Green 
function G a (n 1 z) as 



A a (k,lo) — — Im [Tr G a (n, u> + i5)} /it. 



(4) 



Then the renormalized density of states p a (^) is given by 
Pit(w) = 2N~ 1 A <7 (ek, u>) where the summation runs 
over the reduced Brillouin zone with N/2 points. In the 
paramagnetic state, the reduced number of k is compen- 
sated by the trace over Gk a (z) to give the same p a (^>) 
as derived by use of the original Brillouin zone. Even in 
the Neel state, p a (u) does not depend on a because of 
summation over sublattices. 

Let us first discuss the magnitude of the staggered mo- 
ment given by 2(Sq) = (S%) — (S#), where we choose 
the positive polarization for the A sublattice. Figure 
Q] shows the temperature dependence of the staggered 
moment, which should vanish at the Neel temperature 
Tn- On the other hand, the staggered magnetic suscep- 
tibility xq m th° paramagnetic state should diverge as 
the temperature is lowered toward Tn. Hence we also 
plot Xq 1 calculated in ref . [3| . It is found that the es- 
timates of Tn by (Sq) and by Xq 1 are consistent with 
each other. However, calculation of (Sq) becomes in- 
creasingly difficult as T is approached to Tn- The results 
shown in FIGfT] are restricted to the temperature range 
where we could obtain reliable values. With J = 0.1, the 
localized spins are almost fully polarized in the ground 
state. As J increases, (Sq) decreases by the Kondo ef- 
fect, which eventually suppresses the antiferromagnetism 
at J = J c ~ 0.27 down to T = 0. 

We next discuss the density of states near T = 0, 
which is derived using the Pade approximation for ana- 
lytic continuation from imaginary Matsubara frequencies 
ie„ to the real ones. Since the Monte Carlo data are ob- 
tained accurately in the imaginary time domain, the Pade 
approximation works well in the CT-QMC method [l3T|. 



Figure [2] shows p a (u>) for different values of J near the 
ground state. Namely, we take the temperature where 
the density of states does not vary much when T de- 
creases further. For example, we take T = 0.01 for 
J = 0.3, but T = 0.002 for J = 0.1. The present method 
has no difficulty to go to such low temperatures with high 
accuracy. In the case of J = 0.3 with the paramagnetic 
ground state, the density of states has a gap caused by 
the Kondo effect. The state is often called the Kondo 
insulator. In the case of J = 0.26, the ground state is 
antiferromagnetic as seen from FIGfTJ The density of 
states in the ordered phase is almost the same as that 
with J = 0.3. Namely the density of states does not de- 
pend much on whether the ground state is paramagnetic 
or antiferromagnetic as far as J is close to the critical 
value J c . 

When J is smaller, sharp peaks develop at both edges 
of the gap. The origin is understood as follows: Putting 
the self-energy of conduction electrons as the staggered 
potential ±h by the Neel order, we obtain 



TrG a (n, z) 
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z + \J k 2 + h 2 z — \J k 2 + h 2 
Then the density of states is given by 
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Pc (Vu 2 -h 2 ), 



(5) 



(6) 



with a square-root divergence at both edges lo = ±h. 
For lo 2 < h 2 , we obtain p a (d) = 0. This form of p a (u) 
roughly explains the peak structure with J = 0.1 in 
FIGlH In the numerical result, the residual Kondo ef- 
fect actually suppresses the divergence in p a (w). 

On the other hand, the density of states for larger J 
does not show a clear threshold. The Gaussian tail in the 
bare density of states p c (w) causes a tiny but finite magni- 
tude of p a (ijj) within the apparent energy gap. Therefore 
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FIG. 1: (color online) Staggered spin polarization 2(Sq) (left 
scale) as a function of temperature for different values of J. 
Also shown is the inverse staggered susceptibility Xq 1 (right 
scale), which goes to zero at Tn 
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FIG. 4: (color online) Single-particle spectrum with J = 0.2 
in (a) paramagnetic phase at T — 0.035 and (b) antiferromag- 
netic phase at T = 0.010. 



FIG. 2: (color online) J dependence of the density of states 
for the conduction electron and with T ~ 0. 
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FIG. 3: (color online) J dependence of the energy gap A c /2. 
For comparison, also plotted are J{S Z ), J (a*), S z and si are 
the localized and conduction spins at a local site, respectively. 
V 2 /D represents the effective hybridization to be explained 
later, and 7k is the Kondo temperature defined by Eq.©- 



we introduce the characteristic value A c of the energy 
gap as giving the half-maximum value of the peak in the 
density of states. Figure [3] shows A c /2 as a function of J. 
The error bars have been estimated from 5 bins of data. 
It is clear that A c changes continuously at J = J c . This 
shows that both Kondo effect and the staggered internal 
field are contributing to A c . For J < 0.1, on the other 
hand, A c /2 is almost proportional to J. This behavior 
shows that the gap is mainly determined by the staggered 
field J(S Z ) as shown also in FIG. [3]. 

The details of the itinerancy are seen in the single- 
partcle spectral function A CT (/t, ui). Figure H] shows the 
spectrum in the case J = 0.2 where the Kondo effect is 
significant. In both paramagnetic and ordered phases, 
the spectrum consists of four bands in the reduced Bril- 
louin zone. In the paramagnetic phase shown in the 
left panel, the new bands are ascribed to "hybridization 
bands" caused by the Kondo effect [lj|. We note that 
there is no real hybridization between the /-electron and 
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FIG. 5: (color online) A a (n, uS) with ferromagnetic interaction 
J = —0.2 in (a) paramagnetic phase at T = 0.035 and (b) 
antiferromagnetic phase at T = 0.010. 



the conduction electron because the /-electron does not 
have the charge degrees of freedom in the KLM. In the 
Brillouin zone of the paramagnetic state, the energy gap 
is indirect from the zone boundary to the zone center, 
both of which correspond to k — > ±oo in FIG|H In the 
right panel of FIGrjJb), these hybridization bands are 
clearly seen even in the antiferromagnetic phase. Hence 
we classify the antiferromagnetism in this regime as itin- 
erant. 

Let us now present the spectrum in the ferromagnetic 
KLM for comparison, where each site forms S = 1 state 
with antiferromagnetic intersite interaction. The Kondo 
effect is absent in the case of J < 0. We can apply the 
CT-QMC method also to ferromagnetic J as noted in ref. 
17]. Figure [S] shows the spectrum with the same T as in 
FIG|H In the paramagnetic state shown in FIGG^a), the 
spectrum is almost the same as the original conduction 
band. On the other hand, in the antiferromagnetic phase 
shown in FIGOJb), the spectrum shows the clear gap 
structure. These behaviors are explained well by Eq.Q 
with h = in FIGlUJa), while h ^ in (b). In the case of 
ferromagnetic J, the number of bands is two in contrast 
to the case with J > with four bands. 

Thus, the itinerant or localized behavior can be distin- 
guished by the number of energy bands in the reduced 
Brillouin zone. Moreover the location of the energy gap 
is at the zone boundary in the localized case, but at the 
zone center in the itinerant case. Note that k — ±oo 
correspond to the lower and upper edges of the conduc- 
tion band in infinite dimensions. Both edges come to the 
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center in the reduced Brillouin zone. 

We note another characteristic feature in FIG[4|b) that 
the degeneracy at K — remains in the antifcrromagnetic 
phase. A toy model helps to understand the origin of the 
degeneracy. Let us consider a non-interacting periodic 
Anderson model under staggered field as 



HpAM = *^2 £ kcl a C k a + V/Xfi,rC i(T + h.C. 



HiSf), 



(7) 



where fia(fi a ) is the annihilation (creation) operator of 
an /-electron at the i site with the hybridization V. The 
staggered fields represent the molecular field associated 
with the antiferromagnetism. Then we choose hi = ±h 
and Hi = ±H, where the A (B) sublattice has the posi- 
tive (negative) sign. 

There are four branches associated with c and f degrees 
of freedom for the electrons, as well as presence of A and 
B sublattices. At k = 0, we obtain the energies 



l r 



E(k = 0)=±- (h + H)± y/(h - H) 2 + W 2 



(8) 



If the relation h = —H holds, E(k = 0) has only two 
distinct values both of which are doubly degenerate. We 
interpret the degeneracy in FIG|4jb) as caused by the 
relation h = —H for the molecular field. We call such sit- 
uation "quasi-local compensation" . Note that the mag- 
nitudes of the polarization for the conduction spin J{s*) 
and localized spin J(S Z ) are much different as shown in 
FIG El Hence the actual energy level is determined not 
by the local internal field, but by a long-range field in- 
volving remote conduction electrons. The compensation 
is reminiscent of the spatially extended Kondo singlet 
with small J. 

Let us identify the energy scale in the Neel state from 
the self energy. For sublattice a = ±1 and spin a = ±1, 
Y> aa {z) is expanded as 



V 2 ( 1 

^aa(z) = aah H hO — 

z \ z z 



(9) 



where V is the effective hybridization. The coefficient V 2 
of l/z corresponds to the jump at r = in the imaginary 
time domain. We extract numerically the coefficient V 2 
of l/ie„ in the self energy. Figure [3] shows the result for 
V 2 /D (= V 2 ) as a function of J. Note that the indirect 
energy gap in the toy model ([7]) is given by V 2 /D. The 
agreement between V 2 / D and A c /2 in the paramagnetic 
phase is excellent. However, this consistency should not 
be taken too seriously because A c depends on the def- 
inition of the gap. We emphasize that V 2 / D shows no 
anomaly across the phase transition to the antifcrromag- 
netic phase. It also shows good proportionality to the 
Kondo temperature as V 2 / D ~ 2.6Tk. 



In the region < J < 0.1, the Neel temperature Tn 
is much larger than Tk as seen from FIGSQ]and[3l The 
electronic state at the transition has a localized character 
since the Kondo effect is negligible at Tn. Namely there 
is no hybridized band, and the energy gap occurs at the 
boundary of the reduced Brilluoin zone. However, we 
have checked that two almost flat bands appear newly 
below Tk even with J = 0.05. Hence, the crossover from 
localized behavior to the itinerant one occurs inside the 
Neel ordered state. 

In summary, we have derived single-particle spectrum 
and the temperature-dependent order parameter in the 
infinite dimensional KLM allowing for the Neel order. 
The high numerical accuracy has made it possible to 
find the quasi-local compensation between the conduc- 
tion and localized spins indicating the persistent ten- 
dency toward the Kondo singlet at each site. The ef- 
fective hybridization energy has no anomaly across the 
quantum phase transition, and scales well with the impu- 
rity Kondo temperature Tk- Hence the quantum transi- 
tion into antiferromagnetism occurs within the itinerant 
regime, and does not involve itinerant-localized transi- 
tion. 
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